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Abstract. We analyze the space <S^ ^, (7") of bivariato functions that are 
piecewise polynomial of bi-degree ^ (m, m') and of smoothness r along the in- 
terior edges of a planar T-mesh T. We give new combinatorial lower and upper 
bounds for the dimension of this space by exploiting homological techniques. 
We relate this dimension to the weight of the maximal interior segments of 
the T-mesh, defined for an ordering of these maximal interior segments. We 
show that the lower and upper bounds coincide, for high enough degrees or for 
hierarchical T-meshes which are enough regular. We give a rule of subdivision 
to construct hierarchical T-meshes for which these lower and upper bounds co- 
incide. Finally, we illustrate these results by analyzing spline spaces of small 
degrees and smoothness. 



Introduction 

Standard parametrisations of surfaces in Computer Aided Geometric Design are 
based on tensor product B-spline functions, defined from a grid of nodes over a 
rectangular domain [llj . These representations are easy to control but their re- 
finement has some drawback. Inserting a node in one direction of the parameter 
domain implies the insertion of several control points in the other directions. If 
for instance, regions along the diagonal of the parameter domain should be refined, 
this would create a fine grid in some regions where it is not needed. To avoid this 
problem, while extending the standard tensor product representation of CAGD, 
spline functions associated to a subdivision with T-junctions instead of a grid, have 
been studied. Such a T-mesh is a partition of a domain n into axis-aligned boxes, 
called the cells of the T-mesh. 

The first type of T-splines introduced in [211 ED] , are defined by blending func- 
tions which are products of univariate B-spline basis functions associated to some 
nodes of the subdivision. They are piecewise polynomial functions, but the pieces 
where these functions are polynomial do not match with the cells of the T-subdivision, 
Moreover, there is no proof that these piecewise polynomial functions are linearly 
independent. Indeed, [3] shows that in some cases, these blending T-spline func- 
tions are not linearly independent. Another issue related to this construction is 
that there is no characterization of the vector space spanned by these functions. 
For this reason, the partition of unity property which is useful in CAGD is not avail- 
able directly in this space. The spline functions have to be replaced by piecewise 
rational functions, so that these piecewise rational functions sum to 1. However, 
this construction complexifies the practical use of such T-splines. 
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Being able to describe a basis of the vector space of piecewise polynomials of a 
given smoothness on a T-mesh is an important but non-trivial issue. It yields a 
construction of piecewise polynomial functions on the T-subdivision which form a 
partition of unity so that the use of piecewise rational functions is not required. It 
has also a direct impact in approximation problems such as surface reconstruction 
[5] or isogeometric analysis |T4| , where controlling the space of functions used to ap- 
proximate a solution is critical. In CAGD, it also provides more degrees of freedom 
to control a shape. This explains why further works have been developed to under- 
stand better the space of piecewise polynomial functions with given smoothness on 
a T-subdivision. 

To tackle these issues special families of splines on T-meshes have been stud- 
ied. In [6], [8], these splines are piecewise polynomial functions on a hierarchical 
T-subdivision. They are called PHT-splines (Polynomial Hierarchical T-splines). 
Dimension formulae of the spline space on such a subdivision have been proposed 
when the degree is high enough compared to the smoothness [6], |13| . |16| and in 
some cases for biquadratic piecewise polynomial functions [7] . The construction 
of a basis is described for bicubic spline spaces in terms of the coefficients of the 
polynomials in the Bernstein basis attached to a cell. When a cell is subdivided 
into 4 subcells, the Bernstein coefhcients of the basis functions of the old level are 
modified and new linearly independent functions are introduced, using Bernstein 
bases on the cells at the new level. 

In this paper, we analyse the dimension of the space S'^ ^, (T) of bivariate func- 
tions that are piecewise polynomial of bi-degree ^ (m, m') of smoothness r along the 
interior edges of a general planar T-mesh T, where r is a smoothness distribution 
on T- 

As we will see, computing this dimension reduces to compute the dimension of 
the kernel of a certain linear map (namely the map 82 introduced in Section [2|). 

Thus for a given T-mesh, a given smoothness distribution r and a given bi-degree 
(to, to'), it is possible to compute the dimension of 5,^ „j, (T) by linear algebra tools 
(see eg. a software implementation developed by P. AlfecQ for such computations). 
We would like to avoid a case-by-case treatment and to describe this dimension in 
terms of combinatorial quantities attached to T and easy to evaluate. As shown in 
[17j or [I] , the dimension may also depend on the geometry of the T-mesh and not 
just on its topology. This explains why it is not always possible to provide a purely 
combinatorial formula for the dimension of 5^ m'(7"). 

The main results in this paper are 

• a description of the dimension in terms of a combinatorial part that depends 
only on the topology of the T-mesh and an homological part that takes into 
account the fact that the dimension may depend on the geometry of the 
T-mesh (Theorem [3T]) ; 

• combinatorial lower and upper bounds on the dimension that are easy to 
evaluate (Theorem 13. 7p : 

• sufficient conditions under which the lower and upper bounds coincide so 
that the dimension depends only on the topology of the T-mesh (Theorem 
ESI). 
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We proceed as follows. By extending honiological techniques developed in f2l 
and |19j . we obtain combinatorial lower and upper bounds on the dimension of 
these spline spaces for general T-subdivisions. We relate the upper bound to the 
maximal interior segments and their weights and show that the lower and upper 
bounds coincide for T-meshes which are enough regular. Namely, if a T-mesh is 
(m + l,m' + l)-weighted, the dimension depends directly on the number of faces, 
interior edges and interior points. In particular, we obtain the dimension formula 
for a constant smoothness distribution r = (r, r') with m >2r + l and m' > 2r' + 1, 
providing a new proof of a result also available in [51, [T31, [TO for a hierarchical 
T-mesh. The algebraic approach gives an homological interpretation of the method 
called Smoothing Cofactor-Conformality method in [23]. It allows us to generalize 
the dimension formulae obtained by this technique [1^, [13]. We also give a rule 
of subdivision to construct hierarchical T-meshes for which the lower and upper 
bounds coincide. As a consequence, we can recover the dimension of the space 
of Locally Refined splines described in We do not consider the problem of 
constructing explicit bases for these spline spaces, which will be analyzed separately. 

In the first section, we recall the notations and the polynomial properties which 
are needed in the following. Section [2] describes the chain complex associated to 
the spline space and analyzes its homology. In Section [2J we give lower and upper 
bounds on the dimension of the spline space and analyze cases where these bounds 
are coincide. Section |4] deals with the properties of hierarchical T-meshes, obtained 
by recursive subdivisions of cells. In the last section, we analyse some examples for 
small degree and smoothness. 

1. Planar T-splines 

In the following, we will deal with notions which are of topological and algebraic 
nature. We start by the topological definitions. 

1.1. T-meshes. For any set S C IR^, S is its closure for the usual topology, 5° its 
relative interior, dS its boundary. 
We define a T-mesh T of IR^ as: 

• a finite set T2 of closed axis-aligned rectangles of IR^, 

• a finite set Ti of closed axis-aligned segments included in Ua^T^c)'^: 

• a finite set of points To = UreTidr, 

such that 

• For a E 72, da is the finite union of elements of 7i . 

• For (7. (t' G 72 with a a\ a H a' = da D da' is the finite union of elements 
ofTiUTo. 

• For r, t' G Ti with t ^ t' , t D t' ^ dr n dr' C %. 

We denote hy Vl — VJ^^-j-^a C IR'^ and call it the domain of the T-mesh T- 

The elements of T2 are called 2- faces or cells and their number is denoted f2- 
The elements of 71 are called 1-faces or edges. An element of 7i is called an 

interior edge if it intersects Vl° . It is called a boundary edge otherwise. The set of 

interior edges is denoted by T° ■ The number of edges in 7i is /i and the number 

of interior edges is f°. 

An edge parallel to the first (resp. second) axis of is called horizontal (resp. 

vertical). Let 7^°''' (resp. T|°'^) be the set of horizontal (resp. vertical) interior 
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edges and (resp. /{') the number of interior horizontal (resp. vertical) edges. 
Then, the number of interior edges is /f = /{* + /i • 

The elements of To are called 0- faces or vertices. A vertex is interior if it is in f2°. 
It is a boundary vertex otherwise. The set of interior vertices is denoted 7^°. We 
denote by /o be the number of vertices of To and by f§ be the number of interior 
vertices. 

A vertex is a crossing vertex if it is an interior vertex and belongs to 4 distinct 
elements of Ti. A vertex is a T- vertex if it is an interior vertex and belongs to 
exactly 3 distinct elements of 7i. Let (resp. Jq ) be the number of crossing 
(resp. T) vertices. A boundary vertex is a vertex in To H dQ. The number of 
boundary vertices is /q. A vertex is a corner vertex if it belongs to dfl and to a 
vertical and a horizontal boundary edge. 

To simplify the definitions and remove redundant edges, we will assume that any 
vertex 7 G To belongs at least to one horizontal edge Th{'y) G Ti and one vertical 
edge t^{j) G Ti- 

We denote by VhiT) — {si, . . . , s;} C IR (resp. Vv{T) = {ti, . . . , tm} C R) the 
set of first (resp. second) coordinates of the points in vertical (horizontal) segments 
G Ti (resp. G Ti)- The elements of Vh{T) (resp. Vy{T)) are called the horizontal 
(resp. vertical) nodes of the T-mesh T. 

Example 1.1. Let us illustrate the previous definitions on the following T-mesh: 
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In this example, there are /2 = 7 rectangles, /f = 9 interior edges such that = 4 
are horizontal and /i = 5 are vertical. There are /q = 3 interior points 71,72,73; 
71,73 are T- vertices and 72 is a crossing vertex. There are /q = 15 boundary 
vertices and 12 corner vertices. 

The horizontal nodes are Vh{T) = {0, . . . , 5} and the vertical nodes are Vv{T) — 
{0,...,4}. 

For our analysis of spline spaces on T-meshes, we assume the following: 
Assumption: The domain D, is simply connected and n° is connected. 

This implies that ft has one connected component with no "hole" and that the 
number of boundary edges through a boundary vertex is 2. 

1.2. T-splines. We are going now to define the spaces of piecewise polynomial 
functions on a T-mesh, with bounded degrees and given smoothness. An element 
in such a space is called a spline function. 

Definition 1.2. A smoothness distribution on a T-mesh T is a pair of maps (r/j, r^) 
from (i//,(r) X VviT)) to N X N. 
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For convenience, we will define the smoothness map r on Ti as follows: for any 
T £ Ti (resp. T e Ti), r(r) = Yh{s) (resp. r(T) = r„(t)) where s £ Vh{T) (resp. 
t £ (T)) is the first (resp. second) coordinate of any point of r. We will also define 
the horizontal and vertical smoothness on To as follows: for any 7 = (s,t) £ To, 
rh(7) = and v^{'^) = Ty{t). 

For r, r' £ N, we say that r is a constant smoothness distribution equal to (r, r') 
if Vs £ MT),r{s) ^r,yt£ y^{r),r{t) = r'. 

Let R = R[s,i] be the polynomial ring with coefRcients in R. For m,m' £ N, 
we denote by Rm,m' = R{m.m') the vector of polynomials in R of degree < m in s 
and < m' in t. An element of Rm,m' is of bi-degree ^ (m, m'). 

The goal of this paper is to analyse the dimension of the space of splines of 
bi-degree < (m, m') and of smoothness r on a T-mesh T, that we define now. 

Definition 1.3. Let T be a T-mesh and r a node smoothness distribution. We 
denote by S'^ ^, (T) the vector space of functions which are polynomials in Rm.m' 
on each cell a £ T2 and of class r(r) in s (resp. in t) at any point of t Oil" if r is 
a vertical (resp. horizontal) interior edge. 

We will say that / £ ^, (T) is of (continuity) class C on T. We notice that 
the boundary edges and their smoothness arc not involved in the characterization 
of a spline function. 

Example 1.4. We consider again the T-mesh of Example ll.il If we take the node 
smoothness distribution r;i(l) = 1, rh(2) = 0, r/i(3) = r;i(4) = r;i(5) = 1, and r„ 
constant equal to 1, then iS| 3(T~) is the vector space of bicubic piecewise polynomial 
functions on T which are in s for s < 2 and s > 2, continuous for s — 2 and 
in t in Q". 

1.3. Polynomial properties. We recall here basic results on the dimension of the 
vector spaces involved in the analysis of 5^ ^, (T) : 

For any r G Ti, let Z,- G i? be a non-zero polynomial (of degree 1) defining the 
line supporting the edge t. Let A'"(t) = ^t'^'^^- We denote by 3'"(t) — (A''(t)) 
the ideal generated by the polynomial A""(t) £ R and by 3'^ ^'(t) = J''(t) ni?m,m' 
its part of bi-degree < {m,m'). Notice that 3''(t) defines the line supporting the 
egde r with multiplicity (r(r) -1-1). By definition, two horizontal (resp. vertical) 
edges Ti,T2 which share a point define the same ideal CT''(ri) = 3'^{t2). We define 
the bi-degree S for any egde r G Ti as follows: 

• S{t) — (r(T) + 1,0) if T is vertical, 

• (5(t) = (0,r(r) -I- 1) if r is horizontal. 

Let 3^(7) — J'^iTy) + 3'^{Th) = {A'^ (Ty) , A'^ {Th j) where Ty,Th £ Ti are vertical 
and horizontal edges such that Ty HTh — {7}. The ideal 3"" (7) defines the point 7 
with multiplicity (r,,(7) + 1) x (r„(7) + 1). We denote by 3^ ^,(7) = a^^„,(T„) + 
m'i'^h)- Notice that these definitions arc independent of the choice of the vertical 
edge Ty and horizontal edge Th which contain 7. The bi-degree of a vertex 7 S To 
is(5(7) = (rft(7)-f l,r,(7)-|-l). 

Here are the basic dimension relations that we will use to analyse the spline 
functions on a T-mesh. 

Lemma 1.5. 

• dim i?m,m' = (to + 1) X (to' 4- 1). 
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• dim(i?„,„//3^_„,(r)) = | 



(to + 1) X (niin(r(T),TO') + 1) if r G Ti 
(niin(r(T), to) + 1) x (to' + 1) if r G Ti 
• dim (i?™,„'/a^ „,(7)) = (min(r/j(7),m) + l) x (min(r,„(7), m') + 1) for all 
1 ^To- 



Proof. To obtain these formulae, we directly check that 

• a basis of Rm,m' is the set of monomials s'^P with 0<i<to, 0<j< to'; 

• a basis of Rm,m' /'S^ m'i'^) the set of monomials s'^P with < z < to and 
< i < min(r(r), to') if r e (rasp. < i < min(r(r), to), < j < to' if 

• a basis of Rm.m' /'^In m'il) the set of monomials s^P with < i < 
min(r;i(7), to), < j < min(ri,(7), to'). 

since the ideal of an edge r e 7i is up to a translation (s''('^'+^) or (f''('^^+^). □ 

An algebraic way to characterise the C""-smoothncss is given by the next lemma: 

Lemma 1.6 ( 2 ). Let r G 7i and let pi,p2 be two polynomials. Their derivatives 
coincide on t up to order r(T) iff pi — P2 G ^''(r). 

In the following, we will need algebraic properties on univariate polynomials. 
We denote by J7 = R[u] the space of univariate polynomials in the variable u with 
coefficients in IR. Let Un denote the space of polynomials of U of degree < n. For 
a polynomial g € U of degree d and an integer n > d, g Un-d is the vector space of 
multiples of g which are of degree < n. For polynomials gi, . . . , gk ^ U respectively 
of degree di, . . . ,dk and an integer n > max^^i ...^^ di, X^iLi 9i Un-di is the vector 
space of sums of multiples of gi of degree < n. 

We will use the apolar product defined on C/„ by: 



(/, 9)n X! ( i" ) 

where / — X]"=o /^'"^ 9 — S"=o5»^' ^ IRiuln - One of the properties that we will 
need is the following [TH], [H], [TU] : 

Lemma 1.7. Let g G Un, d < n and a G R. Then g is orthogonal to {u ~ a)'^Un-d 
for the apolar product iff 

d^g{a) =0,fc = 0,...,n-d. 
Proposition 1.8. Let ai, . . . , a; G R 6e Z distinct points and di, . . . , G N . Then 

dim j ^^(u — ai)''''Un-di ] = min(n + 1, n — di + 1). 
\i=i / 1=1 

Proof. In order to compute the dimension of V ~ o-iYUn-d CI Un, 

compute the dimension of the orthogonal in [/„ of V for the apolar product. 
Let g G Un he an element of the orthogonal V-^ of V. By lemma fTTTl d^g{ai) = 
0,k = 0, . . . ,n — di,i = 1, . . . ,1. In other words, g is divisible by {u — ai)"~'''+^ for 
i = 1, . . . , / . As the points ai are distinct, g is divisible by 

i=l 

Conversely, any multiple of 11 of degree ^ n is in V-^. Thus = (11) n Un- 
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This vector space V-^ of multiples of 11 in degree n is of dimension max(0,n + 
1 — deg(n)), so V is of dimension 

I 

n + 1 — max(0, n + 1 — deg(n)) = min(ri + 1, deg(n)) = min(n + 1, ^ n — di + 1). 

1=1 

□ 



We are going to use an equivalent formulation of this result: 



(1.1) dim Un/ Y,(^ ~ ai)iUn-d^ = (n + 1 - ^ n - + 1)+ 

\ 1=1 / i=l 

where for any a £ Z, a+ = max(0, a). A similar result is proved in |16) [Lemma 
2] when all di are equal to d, by analyzing the coefficient matrix of generators of 



1.4. Maximal interior segments. In order to simplify the analysis of S^-^ m'(T), 
we introduce the following definitions: 

For any interior edge r e ^ we define p(r) as the maximal segment made of 
edges e Ti of the same direction as t, which contains r and such that their union is 
connected. We say that the maximal segment p(r) is interior if it does not intersect 
the boundary of fJ. 

As all the edges belonging to a maximal segment p have the same supporting 
line, we can define A'"(p) — A'"(t) for any edge r belonging io p = p{t). 

The set of all maximal interior segments is denoted by Mis(T). The set of 
horizontal (resp. vertical) maximal interior segments of T is denoted by MiSft(7~) 
(resp. MIS„(T)). 

The degree of p G Mis(T) is by definition 5{p) — 5{t) for any t d p. 

For each interior vertex 7 G 7^°, which is the intersection of an horizontal edge 
Th & Ti with a vertical edge Ty e T", let phij) (resp. Pv{l)) is the corresponding 
horizontal (resp. vertical) maximal segment. We denote by ^\{'^) (resp. AJ^(7)) 
the equations of the corresponding supporting lines to the power r^(7) + 1 (resp. 
r.(7) + l)- 

Notice that the intersection of two distinct maximal interior segments is either 
a T-vertex or a crossing vertex. 

We say that p S Mis(T) is blocking p' G Mis(T) if one of the end points of p' is 
in the interior of p. 

Example 1.9. In the figure below, the maximal interior edges are indicated by 
plain segments. 
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In this example, pi is blocking p4 and p2 is blocking p^. 

2. Topological chain complexes 

In this section, we describe the tools from algebraic topology, that we will use. 
For more details, see eg. [H], [T^ . 



2.1. Definitions. We consider the following complexes: 





,(T°) 



m.m' 



,(T°) 







©rerf ['''l-'m,™,' (''') 

i 

©reTf [''']-^™,"i'/-J'm,m' (''') 
i 







; 
; 

i 

©^Sro°W^m,™'/3'L.,n' (7) 





The different vector spaces of these complexes are obtained as the components 
in bi-degree ^ (m,m') of i?-modules generated by (formal) independent elements 
[a], [t], [7] indexed respectively by the faces, the interior edges and interior points 
of T. An oriented edge r e 7i is represented as: [r] = [ab] where a, 6 g To are the 
end points. The opposite edge is represented by [ba]. By convention, \ba\ = — [ab]. 
The maps of the complex d\m,m'{T°) are defined as follows: 

• for each face cr e 72 with its counter-clockwise boundary formed by edges 
n = aia2, .. .,Ti= flifli, 92(ct) = [ri] © ■ • • ® [t;] = [0102] ® • • • ® [ajoi]; 

• for each interior edge t € from 71 to 72 G To, c^i(H) = [72] — [71] where 
[7] = if 7 ^ V; 

• for each interior point 7 G Tq, da{['y]) — 0. 

By construction, we have di o di+i = for i = 0,1. The maps of the complex 
'^m m'(7~°) obtained from those of yim.m'{T°) by restriction, those of the com- 
plex S^^, (T~°) are obtained by taking the quotient by the corresponding vector 
spaces of 3mm/(T~°). The corresponding differentials of the complex are denoted 

For each column of the diagram, the vertical maps are respectively the inclusion 
map and the quotient map. 

The complex ^Hm,m'(T~°) is also known as the chain complex of T relative to its 
boundary dT- 
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Example 2.1. We consider the foUowing subdivision T of a rectangle fi: 

H il l 1 1 1. 1 ii.mmjj,LL | ,im.LLi i ijimi i i.i L | iim.ii 



7i 



0-2 



0-3 



/33 



02 



We have 

• 92([ai]) = [7i/3i] + [/337i],52([a2]) 
[/327i], 

• 5i([/3i7i]) = [7i],5i([/327i]; 

• 9o([7i]) =0. 

This defines the following complex: 

i'T) • 0i^i[o-i]i?m,m' ^ ®,j^i[ft7l]^m,m' [7l]^m,r^ 

The matrices of these maps in the canonical (monomial) bases are 



71/32] + [/3i7i],92([a3]) = bi/Sa] + 
7i]!5i([/337i] = [71], 

^ 



[^2] = 




[9i 



III 



where / is the identity matrix of size {m + 1) x (m' + 1) (ie. the dimension of 

Rm^m' ) ■ 

Let US consider the case where 71 = (0, 0), (m, m') = (2, 2) and r is the constant 



distribution (1, 1) on T- The matrices of the complex 62'2('^) ^'^^ 
/ -[Hi] [Hi] \ 

-[n2] [n2] Adi] = { [Pi] [P2] [Pa] ) 

V [Ha] -[Ha] / 
where [11;] (resp. [Pi\) are the matrices of the projections 



Hi = Ha : i?2,2 ^ ^2.2/(52) 
p ^ p mod 



n2 : i?2,2 ^ P2,2/(i') 

p 1^ p mod 



Pi = Pa : i?2,2/(52) ^ R2,2/{s^t^) P2 : i?2,2/(t') ^ P2,2/(S^^') 

p mod i-> p mod (s^,t^) p mod 1— > p mod (s^,t^). 

The matrices 11; are of size 12 x 16 and the matrices P; are of size 9 x 12. 

2.2. Their homology. In this section, we analyse the homology of the different 
complexes. The homology of a chain complex of a triangulation of a (planar) 
domain is well-known [22] [Chap. 4], [12j [Chap. 2]. Since, it is not explicit in the 
litterature, we give in the appendix a simple proof of the exactness of 9^m,m'(T°). 
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2.2.1. The 0-homology. Wc start by analysing the homology on the vertices. 
Lemma 2.2. i/o(9^m.™' (T°)) = -ffo(6™,™'(r°)) = 0. 

Proof. By Proposition ID . 1 1 in the appendix, we have i^o(^m,m' (7"°)) = 0. Taking 
the quotient by ^, (ri), we still get that d\ is surjective so that i?o(©m m' C^")) — 
0. ' ' □ 

Let us describe in more details 

i?oP^;,™.(r°)) = b]T{j)/d,{ [r]3-(r)). 

We consider the free i?-module generated by the half- edge elements [j\t], for all 
interior edges t ^ and all vertices -f ^ t. By convention [j\t] = if 7 G dfl. 

For 7 e 7^°, let Eh{-f) (resp. -£^(7)) be the set of horizontal (resp. vertical) 
interior edges that contain 7 and let E{'y) = Ehij) U £^1,(7). We consider first the 
map 

(7) 

re£;(7) 

[7|r] ^ [7]A--(r) 

By definition of 3J„ m'(7)' t'^i^ surjective. Its kernel is denoted m'il)- 

Let Ph{j) (resp. -Pu(7)) be the set of pairs (r, r') of distinct horizontal (resp. 
vertical) interior edges which contain 7 (with (t, r') identify to (r',r)). We denote 
by P(7) = Phil) U -Pu(7)- If 7 is a T-junction, one of the two sets is empty and 
the other is a singleton containing one pair. If 7 is a crossing vertex, each set is a 
singleton. 

The following proposition gives an explicit description of the kernel „j, (7) , 
that we will exploit hereafter. 

Proposition 2.3. 

Ki,^'il) = ([7k] - [7k'])-R(m,m')-'5M 

(r,r')eP(7) 

+ ([tI^] ^'"(^') - [tI^'] A'-(r))i?(,„„,_i,„,,_,,_i) 

Te-Bfe(7),7-'6-B„(7) 

Proof. Let us suppose first that 7 is a crossing vertex. We denote by ti , T2 the 
horizontal edges, t^^t^ the vertical edges at 7. The matrix of the map ip-^ in the 
basis [jln] is 

[(^,1 = ( A A A' A' ) 

where A = A'"(ri) = A'"(rJ, A' = A^ir^) = A'"(r4). Since A and A' have 
no common factor, the kernel of the matrix [(py] is generated by the elements 
[7ki] ~ [7^2], [71^2] A' - [7|t3]A, [7|r3] — [7IT4], which give the description of 
■^m m' (7) bi-degree < (m, m'). A similar proof applies when there is no horizontal 
or vertical pair of distinct edges at 7. This proves the result. □ 

We use the maps (ipj) ^(zj-^^ to define 

<^ : [7k] -R(m.m')-5(-) ^ [7] •^m.m' (7)1 

76X0° reB^ 7ero° 
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SO that we have the following exact sequence: 

^ K.,m'{l) [7k] Ri^,n')-8(r) ^ [l] ^^^^-{l) ^ 

Using this exact sequence, we can now identify ^^^q-o [7] 3^^ (7) with the quotient 

[7|T]i?(m,m')-5(T)/ XI ■*^™.'"'(^)- 

The next proposition uses this identification and Proposition 12.31 to describe 
more explicitly -ffo(3™,,„' (7"°)): 

Proposition 2.4. We have 

HoiT^,m'(.r'')) =00 b\r]Rim,rn')-Sir) 

([7k] - m,m')— (5(r) 

^(r,r')eP(7) 
+ ([^1^] ~ [7'k])^(m,m')-'5(r) 

+ ([7|T]A"^(r') - [7|r']A'-(r))i?(„,„,)_,(^) 

TG-Bfe(7),i-'e-B„(7) 

Proof. The application 

5i: 0W3^,™'(t) ^ 0[7]3^,™'(7) 
rer° 7^7;° 

lift to an application: 

dl ■■ [T]R{.m,rn')~5{T) ^00 [lh]R{rn,rn')-5iT) 

T [7|r] - [7'|t] 

so that the image of di lift in 0^e7^o 0reE(7)[7k]^(m,m')-5(T) to 

im^i = X! ([^1'^] ^ [^'I'^])^("."')-<5(t)- 

Consequently, 

^^o(©:;,™'(r°)) - [7|r]i?(„,™,)_,(,) / [imai+ X J?^,„,(7) 

76ro°rGi3(7) / \ 7ero° 

which yields the desired description of i?o(©m m'C^))- '-' 

In the next proposition, we simplify further the description of 77o(3m m'0'°))- 
Proposition 2.5. 

^o(3m,„i'(T°)) = ®pGMIs(r) [P] -R(m,m')-5(p) 

/ (E7Gr„°([p.(7)] ^lil) ~ [Phh)] Ar(7))i?(„,„,)_,(^)) . 
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Proof. Let B = ®^^r,r ®re£;(7)[7k] ^(m,m')-5(r), K = imai + 'E-yeT- K^i.rn'il) 

and 

K' = (E(r,r')eP(7)([7k] - [7k'])-R(m,m')-'5(r) 

+ Er=(7,7')eri°([7k] - [7'k])^(m,m')-5(r)) ■ 

As K' C K C B, we have B/J-iT = {B/K')/{K/K'). Taking the quotient by K' 
means, that we identify the horizontal (resp. vertical) edges which share a vertex. 
Thus all horizontal (resp. vertical) edges which are contained in a maximal segment 
p of T are identify to a single element, that we denote [p]. As [7|r] = if 7 £ dil, 
we also have [p] = if the maximal segment p intersects the boundary dfl. This 
yields the desired description of i?o(^m m'C^))- '-' 

Definition 2.6. Let K,,,Ar) = dim Ho{3Z'^,{T")). 

2.2.2. The 1-homology. We consider now the homology on the edges. We use the 
property that -ffi($Hm,m' (T°)) = (see Proposition ID . 21 in the appendix). 

Proposition 2.7. i?i(6^_„,(r°)) = i?o(3^,™,(r°)). 

Proof. As i?o(9^m,m'(T°)) = and i?i(5K™,„'(r°)) = 0, we deduce from the long 
exact sequence (see Appendix IbI) 

• • • ^ i?i(9i™,,„'(r°)) ^ i/i (©;,,„, (r°)) ^ iio(3™,™'(r°)) ^ i/o(%«,™'(r°)) ^ • • • 
that Fi(6j;,,,„,(r°)) ^ i/o(a^„,„,(r°)). □ 

2.2.3. The 2-homology. Finally, the homology on the 2-faces will give us informa- 
tion on the spline space ^, (T) . 

We have the following result (proved in Proposition ID. 31 in the appendix): 

Proposition 2.8. ff2(«m,m'(T°)) = Rm,m'. 

The following proposition relates the spline space S^^ ^, (T) to an homology 
module. 

Proposition 2.9. i72(6^^„,(r°)) - ker^a - Sl,^^,{T). 

Proof. An element of H2 (©^ (7"°)) = ker 82 is a collection of polynomials {pa)a£T2 
with G Rm,m' and Po- = p^-' mod CJ'5^(t) if cr and ct' share the (internal) edge r. 
By Lemma 11.61 this implies that the piecewise polynomial function which is on 
(T and Pa' on a' is of class C across the edge r. As this is true for all interior edges, 
ij>a)a£T2 G kerSj is a piecewise polynomial function of Rm,m' which is of class C, 
that is an element of S^ ^,{T). □ 

3. Lower and upper bounds on the dimension 
In this section, are the main results on the dimension of the spline space iS^'jj „i'(7'). 
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Theorem 3.1. Let T be a T-mesh and let r be a smoothness distribution on T. 
Then 

(3.1) dim5^,,„,(r) = ^(m+l)(m' + l) 

aeT-2 

- (™+l)(r(r) + l)- K + l)(r(r) + l) 

+ 5](r,(7) + l)(r.(7) + l) 

where h^„,^^,iT) = dimi/o(3^,,„,,(r°)). 
Proof. The complex 

m,,m' ~ m.m' m.m' rn,,m' m,m' 

(7)^0 

induces the fohowing relations 

^ dim{H2{e^ra,n'iTl)) - A\u,{H,{&^^^, {T''))) + dim(i7o (6^,™- (T"))) 
As i?2(6^„,„,(r°)) = 5,'-„^„,(r), i/o(©5;.,„-(r°)) = O and H^i&^^m'iT'^)) = 
^o(3™,„'(T°)), we deduce that 

dini5;„,„,(r) = dini(©„er2Hi?™^™') -dim(®,eri°H^m.™'/^™,™'M) 
+ dim(©^gro°[7]i?m,m'/3j;,„.(7)) + dim(i7o(a:;;':„,(r))) 
which yields the dimension formula (j3.ip using Lemma ll.51 □ 

As an immediate corollary of this theorem and of Proposition [2?5l we deduce the 
following result: 

Corollary 3.2. // the T-mesh T has no maximal interior segments then h'^ ^, (T) = 
0. 

In the case of a constant smoothness distribution, Theorem 13.11 is written as 
follows: 

Theorem 3.3. Let T be a T-mesh and let r = (r, r') be a constant smoothness 
distribution on T . Then 

(3.2) dim5;;_^,(r) = (m + l)(m' + l)/2 

- ((m + l)(r' + 1)/^ + (m' + l)(r + 1)/^) 
+ (r + l)(r' + l)/o 
+ Kn.m'iT)- 

where 

• /2 is the number of 2- faces G 72, 

• fi (resp. fi) is the number of horizontal (resp. vertical) interior edges 

• /o is the number of interior vertices £ 7^°, 

• /i^,,„,(r)-dimi?o(a^,„.(r°)). 

We are going now to bound h^^ ^, (T) for general T- meshes. 
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Definition 3.4. Let l be an ordering of Mis(T) that is a map from Mis(T) to N. 
For p e Mis(T), let r, (/9) be the set of vertices 7 of p which are not on a maximal 
interior segment p' e Mis(7~) with i(p') > The number of elements of rt(p) is 
denoted At(p). 

We define now the weight of a maximal interior segment. 
Definition 3.5. For p e Mis(T), let 

• ^^^ip) ^ E7Gr, ("^ - ^vh)) if P e Mis,,(r). 

• uj.ip) = E^er.K - T^hh)) if P e Mis.„(r). 
We called it the weight of p. 

As in the usual spline terminology, for an interior point 7 £ Tq, we call 7 
m — r/i(7) (resp. vn! — r/i(7) ) the horizontal (resp. vertical) multiplicity of 7. 

If p is horizontal (resp. vertical), the weight of p is the sum of the vertical (resp. 
horizontal) multiplicities of the vertices 7 G V ^{p\ 

Notice that if r = (r, r') is a constant smoothness distribution on T, then ^^(p) = 
(to - r)A(p) for p £ MiSh{T) (resp. ui^ip) = (m' - r')X{p) for p e MiSt,(T)). 

Example 3.6. We consider the T-mesh of Example 11.91 with {m,m') ~ (2,2), 
the constant smoothness distribution r — (1,1) and the ordering of the maximal 
interior segments L{pi) = i for i — I, ... ,4. Then we have 



• 


Wt(pi) 


= (2- 


1) 


X 2 


= 2, 


• 


Wt(p2) 


= (2- 


1) 


X 2 


= 2, 


• 




= (2- 


1) 


X 3 


= 3, 


• 


W,,(p4) 


= (2- 


1) 


X 3 


= 3, 



since the multiplicity of a vertex is 2 — 1 = 1 and the interior points of pi , p2 are 
not in T,,{pi) or T,{p2). 

In the following, we will drop the index t to simplify the notations, assuming 
that the ordering t is fixed. Then we have the following theorem: 

Theorem 3.7. Let T be a T-mesh and let r be a smoothness distribution on T- 
Then 

Q < hl,^^,{T) < (™ + l-^(p))+ X (TO'-r(p)) 

peMiSfe(r) 

+ i-(/5)) X (to' + 1 - u{p)))^ . 

pGMIS„(T) 

Proof. Let pi, . . . , pi be the maximal interior segments of 7". By Proposition 12.51 
Knm'(T) is dimension of the quotient in bi-degree < {m,m') of the module 
M := ©■^i[pi] R by the module K generated by the following relations: for each 
vertex 7 G Tq which is on a maximal interior segment, 

• A'^{pj)[p.] — A'^{pi)[pj] if 7 is the intersection of the maximal interior seg- 
ments Pi and pj, 

• A'^(p)[p.] if 7 is the intersection of the maximal interior segment pi with 
another maximal segment p intersecting dfl. 

To compute the dimension of M/K in bi-degree < {m,m'), we use a graduation 
on M given by the indices of the segments. For r :— J2iPi[Pi] £ ^ (with pi G 
R(m.,m')-S{pi))j l6t In(TO) be the element pi„ [pi„] where io is the minimal index such 



ON THE DIMENSION OF SPLINE SPACES ON PLANAR T-MESHES 



15 



that Pi 7^ 0. Wc denote it by In(r) and called it the initial of r. Let In(A') = 
{In(fc) I k e if}. The dimension /i^ m'i'T') is then 

^™,™'(r) = dim{M/K) ^ dim{M/ln{K)). 

Notice that ln{K) contains the multiples in bi-degree ^ (m, m') of 

• A'^{pj)[p.] if 7 is the intersection of two maximal interior segments pi and 
Pj with i > j, 

• A'^{p)[p.] if 7 is the intersection of the maximal interior segment pi with a 
maximal segment p intersecting Q. 

Let Li be the vector space spanned by these initials in bi-degree ^ {m,m'), which 
are multiples of [pi]. By definition, for each 7 G r(pi), we have a generator A'^{p)[pi] 
in Li for {7} — pi f] p. By Proposition II. 8[ Li is of dimension 

• mm{m+l,uj{pij) x (m' - r{pi)) if pi G MIS,,(T), 

• min(m' + l,U}{p,)) x (m - r{pi)) if pi G MlS„(r). 
Thus the dimension of [pi]R(m,m')-S{pi)/Li is 

• (to + 1 - uj{pij)+ X (m' - r{pi)) if p, G MlS,i(r), 

• (to - r{pi)) X (to' + 1 - Lj{p,))+ if p., G MlSt,(T). 
As ln{K) D J2i Li, we have 

h';n,m'{r) = dim(M/In(X)) 

< dim([pj R{rn,m')~S{pi}/(Y^ Li)) = ^ dim ([pJ R[ni.m')-S(pi)l Li) . 

i i 

This gives the announced bound on (T) , using the previous computation of 

dim([pj i?(„^„/)_5(p^)/Lj). □ 

Definition 3.8. The T-mesh T with a smoothness distribution r is (A;, fc')-weighted 
if 

• Vp G Mis,i(T), w(p) > k 

• Vp G MiSu(r), w(p) > k'; 

Theorem 3.9. Let T be a T-mesh and let r be a smoothness distribution on T. If 
T is (to + 1,to' + 1) -weighted, then m'(T) — 0. 

Proof. By definition, Vp G MiSh(T), t^(p) > to + 1 (i.e. (to + 1 — cj(p))+ = 0) and 
Vp G MIS^(T), uj{p) > to' + 1 (i.e. (to' + 1 - uj{p))+ = 0). By Theorem [3T71 we 
directly deduce that h'^^ m'O') — '-' 

Here is a direct corollary which generalizes a result in |16| : 

Corollary 3.10. Suppose that the end points of a maximal interior segment p G 
Mis(7") are in r(p). If for each horizontal (resp. vertical) maximal interior segment 
p G Mis(7"), the sum of the vertical (resp. horizontal) multiplicities of the end points 
and of the vertices of p on a maximal segment connected to the boundary is greater 
than or egual to m \ (resp. m' -\- 1), then h'^ m'iT) = 0- 

Proof. Let p be a maximal interior segment of T. By hypothesis, the end points of p 
are in r(p). As any point 7 G p which is also on a maximal segment connected to the 
boundary is in r(p), the hypothesis implies that oj{p) > to-|- 1 if p is horizontal and 
uj{p) > to' + 1 if p is vertical. We deduce by Theorem 13.71 that h'^ „i'(7") = 0. □ 

Another case where /i^ ^, (T) is known is described in the next proposition. 
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Proposition 3.11. // Vp G MlS/i(T), uj{p) < m + 1 and Vp G MIS«(T), uj{p) < 
m' + 1, then 

(3.3) /i^,™,(r) - (™ + l-^(p))+x(™'-r(p)) 

pGMis,,(r) 

+ {m-r{p)) X {m' + l-u{p)))^. 

peMis„(T) 

Proof. In the case where Vp S MiSh{T), uj{p) < m + 1 and Vp £ MIS„(T), uj{p) < 
m' + 1, Proposition 11.81 imphes that there is no relations in bi-degree ^ {m,m') 
of the monomial multiples of A'^{pj)[p.], A'^{p)[p.] for i — 1, . . . ,l,j < i using the 
same notation as in the proof of Theorem 13.71 This implies that In(if) = ®iLi, 
which shows that 

h'm,m'{r) = dim(M/In(i^)) - ^ dim(i?(„,„,,)_5(p,) [pJ/L,). 

i 

Thus the equality ([331) holds. □ 
As a corollary, we have the following result for constant smoothness distribution: 

Theorem 3.12. Let T be a T-mesh and let r = (r, r') he a constant smoothness 
distribution on T . Then 

(^<K,^^,{T) < Y (™ + l-(™-r-)A(p))+ x(m'-r') 
peMiSh(r) 

+ ^ (m — r) X (m' + 1 — (m' — r')X{p))^ . 

Moreover, equality holds in the following cases: 

• Vp G MlS/i(T), (m - r)A(p) > m + 1 and Vp G MlSi,(T), (m' - r')\{p) > 
m' + 1; 

• Vp G MlS;i(T), (to - r)\{p) < TO + 1 BTid Vp G MIS„(T), (to' - r')A(p) < 
to' + 1. 

4. Hierarchical T-meshes 

We consider now a special family of T-meshes, which can be defined by recursive 
subdivision from an initial rectangular domain $7. Their study is motivated by 
practical applications, where local refinement of tensor-product spline spaces are 
considered eg. in isogeometric analysis |14) . 

Definition 4.1. A hierarchical T-mesh is either the initial axis-aligned rectangle 
ri or obtained from a hierarchical T-mesh by splitting a cell along a vertical or 
horizontal line. 

A hierarchical T-mesh will also be called a T-subdivision. It can be represented 
by a subdivision tree where the nodes are the cells obtained during the subdivision 
and the children of a cell a are the cells obtained by subdividing a. 

In a hierarchical T-mesh T, the maximal interior segments are naturally ordered 
in the way they appear during the subdivision process. This is the ordering l that we 
will consider hereafter. Notice that a maximal interior segment p is transformed by 
a split either into a maximal segment which intersects dVL or into a larger maximal 
interior segment with a larger weight. 

We remark that in a hierarchical T-mesh, if pi is blocking pj then i < j. 
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Example 4.2. Here are an hierarchical T-mesh (case a) and a non-hierarchical 
T-mesh (case b): 



t) — 



(a) 



(b) 



4.1. Dimension formula for hierarchical T-meshes. As a corollary of Theo- 
rem [STTl we deduce the following result, also proved in [6], [16], [T3] : 

Proposition 4.3. Let T be a hierarchical T-mesh and let r = (r, r') be a constant 

smoothness distribution onT. Form > 2r+l andm' > 2r'+l, we have h^^^i^T) = 
0. 



Proof. We order the maximal interior segments in the way they appear during the 
subdivision. If a segment pi E Mis(T) is blocking pj £ Mis(T), we must have i < j- 
This shows that the end points of pi are in T{pi). Thus, X{pi) > 2. 
As TO > 2r -|- 1 , we have 

(to — r)X(pi) > 2{m — r) > m + (to — 2r) > m + 1. 

Thus, (to -I- 1 - (m - r)X{pi))+ = 0. Similarly (m' + 1 - (m' - r')X{p,))+ = holds 
since m' > 2r' + 1. By Theorem gTl we deduce that /i,''„,„'(T) = 0. □ 

Theorem 13.91 leads us to the following construction rule of a T-subdivision T for 
which /j^ ,„,(r)-0. 

Algorithm 4.4 ((fc, /c')-weighted subdivision rule). 
For each 2-face a of a T-mesh to be subdivided, 

(1) Split a with the new edge t; 

(2) // the edge t does not extend an existing segment, extend t (on one side 
and/or the other) so that the maximal segment containing t is either in- 
tersecting dfl or horizontal (resp. vertical) and of weight > k (resp. > k' ). 

If such a rule is applied in the construction of a T-subdivision, 

• either a new maximal interior segment is constructed so that its weight is 
> A; if it is a horizontal maximal interior segment (resp. > k' it it is a 
vertical maximal interior segment), 

• or an existing maximal interior segment is extended and its weight is also 
increased, 

• or a maximal segment intersecting dO. is constructed. 

In all cases, if we start with a {k, fc')-weighted T-mesh, we obtain a new T-mesh, 
which is also (fc, fc')-weighted. 
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By Theorem ESI if A: > m + 1 and fc' > m' + 1 then /i™.„,(T) = and the 
dimension of dimS^^^ ^^^/{T), given by formula (I3.ip . depends only on the number 
of cells, interior segments and interior vertices of T. From this analysis, we deduce 
the dimension formula of the space of Locally Refined splines described in [9] . 



5. Examples 

In this section, we analyse the dimension formula of spline spaces of small bi- 
degree and small constant smoothness distribution r — {r,r') on a T-mesh T. 

5.1. Bilinear C°'° T-splines. We consider first piecewise bilinear polynomials on 
7", which are continuous, that is m = m' = 1 and r ~ r' ^ 0. By Proposition 
14.31 have /ii'i(T) = 0. Using Theorem 13.11 and Lemma [A. II in the appendix, we 
obtain: 

(5.1) dim5°f (T) = 4/2 - 2/° + /° = /+ + ^^ 



5.2. Biquadratic C^-^ T-splines. Let us consider now the set of piecewise bi- 
quadratic functions on a T-mesh T, which are C^. For m = m! = 2 and r = r' = 1, 
Theorem 13.11 and Lemma lA.ll again yield 



(5.2) dim^^^CT) = 9/2 - 6/r + 4/° + h'^^iT) = f+ - ^fj + ^fo + 3 + hl^liT). 



If the T-mesh T is (3, 3)-weighted, then by Theorem 13.91 we have a h2'2{'T) — 0, 
but this is not always the case. 

Example 5.1. Here is an example where h2'2i'T) = 1 by Proposition 13.111 since 
there is one maximal interior segment p with ijj[p) = 2 — 1 + 2 — 1 ^ 2: 



:i4 



to 



S2 



S3 



The dimension of S^^l (r)is9x4-6x5-h4x2-h hl'l{T) = 14-^1 = 15. Notice 
that the dimension is the same without the (horizontal) interior segment. Thus a 
basis of S2 '2 i'T) is the tensor product B-spline basis corresponding to the nodes 
soi soj sq: si, S2, 53, 53, S3 in the horizontal direction and the nodes ^O; ^O; ^i; ^i; ti 
in the vertical direction. 



Example 5.2. Here is another example. We subdivide the T-mesh 7i to obtain 
the second T-mesh J2: 
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0-2 



Doing this, we increase the number of cells by 9 — 1 = 8, the number of interior 
edges by 24 — 4 = 20, the number of interior points by 16 — 4 = 12. The dimension 
of the spline space increases by9x8 — 6x20 + 4x12 



h\'\{T2) — h\'\{Ti). Since there is no maximal interior segment in 7i, by Corollary 
13.21 we have h2\{Ti) = 0. Choosing a proper ordering of the interior segments, 
we deduce by Theorem 13.11 that h\'^{T2) < 1- Suppose that ci — [00,03] x [^c^a] 
and (72 — [01,02] X [61,62]- For uq < ui < U2 < & R, let 7V(m; uq, wi, M2, U3) 
be the B-spline basis function in the variable u of degree 2 associated to the nodes 
Wo, . . . , U3 (see [5]). Then the piecewise polynomial function 

iV(s; oo, oi, 02, 03) X N{t] 60, 61, 62, 63) 
is an element of ^2 '2 (72), with support in cti. It is not in S\'2{Ti), since the function 
is not polynomial on cri. Thus we have dim ^2 '2(72) = dimiS2 2(7i) + 1. 

Notice that 72 is (2, 2)-weighted but not 3-weighted, since any new maximal 
segment intersects two of the other new maximal segments. 

For a general hierarchical T-mesh, we consider a sequence of T- meshes 7o, ■ • ■ 7/ 
where 7o has one cell, 71 = T and such that 71+ 1 refines % by inserting new edges. 
We can assume that at each level i ^ a, new maximal interior segment pi appears 
and that we number the maximal interior segments of T in the order they appear 
during this subdivision. Notice that any maximal interior segment of T extends 
one of the maximal interior segments pi and thus its weight is bigger. Notice also 
that the maximal interior segment pi introduced at level i extends to a maximal 
segment of T, which may intersect the boundary. In this case, it is not involved in 
the dimension upper bound. Then, we have the following corollary: 

Proposition 5.3. Let T be a hierarchical T-mesh. 

(5.3) 9/2 - 6/r + 4/0° < dim Sl;^ (T) < 9/2 - 6/1° + 4/0° + a 

where a is the number of levels of the subdivision where a maximal interior segment 
with no-interior point is introduced. 

Proof. Consider the new maximal interior segment pi of 71 appearing at level i. By 
construction, the end points of pi are not on maximal interior segments of bigger 
index. Thus Lu{pi) > 2. If pi contains an interior vertex then uj{pi) > 3 and 
(3 — CLi(/9i))_|_ = 0. Otherwise {S — Lu{pi))^ = 1. As pi extends to a maximal segment 
Pi which is either interior or intersecting the boundary, we have (3 — aj(pi))+ < 
(3 — aj(pi))+ with the convention that uj{pi) —3 if Pi is intersecting the boundary. 
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By Theorem 13 .71 we have 

I I 
< hllliT) < ^(3 ~ < ^(3 - 

1=1 1=1 
where I is the number of levels in the subdivision, pi is the maximal interior segment 
introduced at level i and pi is its extension in T. By the previous remarks, (3~ 
uj{pi))+ — a is the number of levels of the subdivision where a maximal interior 
segment with no interior point is introduced. Using Theorem 13.11 this proves the 
bound on the dimension of S2'2{'T). □ 

Examples 15.11 and 15 . 2 1 show that the dimension can be given by the upper bound. 
On other other hand, for any (3,3)-weighted hierarchical subdivision, the lower 
bound is reached. This shows that the inequalities (15.31) are optimal for dim52'2 (7~). 

Remark 5.4. In the case of a hierarchical subdivision where some cells of a given 
level are subdivided (as ai in Example 15. 2p into 9 sub-cells which have the same 
length and height, it can be proved that the dimension is in fact: 

dunSl'jiT) = 9/2 - 6/1° + 4/° + a 

where a is the number of isolated subdivided cells (i.e. the cell is subdivided, not 
touching the boundary and the adjacent cells sharing an edge are not subdivided) 
at some level of the subdivision. Indeed, any maximal interior segment subdividing 
a non-isolated cell contains an interior point and is not involved in the upper bound. 
As in Example 15.21 only the isolated cells have a maximal interior segment with 
no interior points. This example also shows that a new C^'^'^"' bi-quadratic basis 
element can be constructed for each isolated cell, proving that the upper bound 
is reached. This gives a dimension formula similar to the one in [7], for a slightly 
different subdivision strategy. 

5.3. Bicubic C^'^ T-splines. For m = m' — 3 and r = r' = 1, that is for piecewise 
bicubic polynomial functions which are C^, Proposition 14.31 vields h^'^{T) = 0. 
Using Theorem 13.11 and Lemma lA.ll in the appendix, we obtain: 

(5.4) dim43(r) = I6/2 - 8/f + 4/° = 4(/+ + f,). 

5.4. Bicubic C^'^ T-splines. For m = m' = 3 and r = r' = 2, by Theorem 13.11 
and Lemma lA.li we have: 

(5.5) dimSl:Ur) = I6/2 - 12/1° + 9/0° + /i^l(r) - /„+ - + 2f, + S + hl^'.iT). 
If 7" is a hierarchical (4, 4)-weighted subdivision, then by Theorem 13.91 we have 

hlim = 0. 
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Appendix A. Combinatorial properties 

We recall some well-known enumeration results for a T-mesh of an axis-aligned 
rectangular domain £7. 

Lemma A.l. 

• h = /o+ + \fl + \f', 1 

. ft = 2/o+ + |/o^ + 2 

• /o° = ./o+ + 

Proof. Each face fi G 72 is a rectangle with 4 corners. If we count these corners for 
all cells in 72, we enumerate 4 times the crossing vertices, 2 times the T- vertices 
which are interior or on the boundary and one time the corner vertices of f7. This 
yields the relation 

4/2=4/+ + 2(/o^ + (/o^-4)) + 4. 

Each interior edge t e 7^° has two end points. Counting these end points for all 
interior edges, we count 4 times the crossing vertices, 3 times the T- vertices which 



22 



BERNARD MOURRAIN 



are interior and one time the T-vertices on the boundary: 

2/r-4/+ + 3/o^ + (/o^-4). 
FinaUy, as an interior vertex is a crossing vertex or a T-vertex, we have 

/o = /o^ + fo- 

□ 

Appendix B. Complexes and homology 

Let us recaU here the basic properties that we wih need on complexes of vector 
spaces. Given a sequence of K- vectors spaces Ai, i = 0,...,/ and hnear maps 
di : Ai ^ Ai+i, we say that we have a complex 

A:Ai^ Ai_i ^---Ai^ Ai_i -^---Ai^Ao 

if imdi C kerdi^i. 

Definition B.l. The i^^ homology Hi{A) of A is ker9j_i/ imdi tor i — 1, ... ,1. 

The complex A is called exact (or an exact sequence) if Hi{A) = (ie. im.di — 
ker9i_i) for i = 1, . . . , /. If the complex is exact and Ai = A^^ = 0, we have 

^(-l)MimA, = 0. 

i=l 

Given complexes A = (Ai)i=o,...,/ ^ ~ (■^i)i=o,...,/, C = (Ci)i=o,...,/ and exact se- 
quences 

Q Ai^ Bi^ Ci^Q 
for i = 0, . . . , we have a long exact sequence [4], [22[p- 182]: 

> H,+i{C) ^ H,{A) ^ H,{B) ^ H,{C) ^ H.^iiA) ^ • • • 

Appendix C. Dual topological complex 

The dual complex of the subdivision T, is such that we have the following 
properties. 

• a face cr G 72 is a vertex of the dual complex T* . 

• an edge of T* is connecting two elements cr, cr' e 72 if they share a common 
(interior) edge t £ Ti ■ Thus it is identified with the edge r of 7" between 
cr, cr'; 

• a face of 7~* corresponds to an elements 7 e . It is either a triangle if 7 
is a T-junction or a quadrangle if 7 is a crossing vertex. 

Notice that the boundary cells of T correspond to boundary vertices of T* . They 
are connected by boundary edges which belong to a single face of T* . 

Appendix D. Topological chain complex 

In this appendix section, we recall the main properties of the topological chain 
complex 

where 

. V7er„°,ao([7]) = o. 

• Vr = [71,72] e 71°, 9i([7i,72]) = [71] - [72] with [7] = iff 7 £ dn. and 
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• Vcr G 72° with its counter-clockwise boundary formed by the edges [71, 72], ... , [7;, 71], 
d2{a) = [71,72] H h [7;, 71] with [7,7'] = iff 7,7' e dCl. 

We assume that fl is simply connected and 0,° is connected. We prove that 

^m.m'{T°) is acyclic on a T-mesh of fi. 

Proposition D.l. -ffo(9^m,m'(T°)) = 0- 

Proof. Let 7 e Tq. There is a sequence of edges tq = 7o7i, n = 7172, 
n = im+i, such that Ti e Ti, 70 ^ Tq and 7^+1 = 7. Then 

di{[To] + ■■■ + [n]) = [71] - [70] + • • • + h+i] - bi] = [7] 

since [70] = and [7;+i] = [7]. Multiplying by any element in Rm,m', we get that 
h]Rm^rn' C im and thus i/o(5^m,m'(T°)) = 0. ' □ 

Proposition D.2. i?i (r°)) 0. 

Proof. Let p = XlTeT°M^'^ ^ keri9i with p^ e Rm,m'- Let us prove that p is in 
the image of 82- For each 7 G Tq" and each edge r which contains 7, we have 
^ErPr = with Et- = 1 if r ends at 7, = — 1 if r starts at 7 and e,- = 
otherwise. 

For any a GT2 and r G T", we define e^.r = 1 if r is oriented counter-clockwise 
on the boundary of cr, e^.r = — 1 if r is oriented counter-clockwise on the boundary 
of a and Eo^t = otherwise. 

For any oriented edge of the dual graph T* from cr' to tr, let us define d*{[a', a]) = 
Ea^rPr- Noticc that 9*([(T, (t']) = Ea'.TPT = —So.tPt = ([^ ') o"])) siucc the orien- 
tation of T on the boundary of a and cr' are opposite. 

Let fJo be the 2-face of T with the lowest left corner for the lexicographic ordering. 
We order the cells cr G according to their distance to erg in this dual graph T* ■ 
We define an element q = Xl<TeT2 9<^['^] where G Rm,m' by induction using this 
order, as follows: 

• qao = 0; 

• For any a, a' G 72 , if cr > cr' and cr and a' share a common edge r, then 
9<T = qa' +dl{[a' ,a]). 

Thus, if [ao, cri], [cri, cr2], . . . , [cr^^i, Gk] is a path of T* connecting erg to cr^ = cr with 
di+i > ai then q„ = J2i=o (^lii'^i^'^i+i])- Let us prove that this definition does not 
depend on the chosen path between ag and cr. 

We first show that for any face 7* of T* attached to a vertex 7, if its counter- 
clockwise boundary is formed by the edges [cr, cr'], [cr', cr"], . . . , [cr'", cr] corresponding 
to the edges t, t', t", . . . of T containing 7, then 

ai*([a, a']) + 9i*([(r', a"] + ■■■ + dl{[a"', a]) = e^,rPr + e^' ,r'Pr' +ea",r"Pr" + • • • = 0. 

By changing the orientation of an edge r, we replace Pr by —pr and e^.r by — e<T,r 
so that the quantity e^.tPt is not changed. Thus we can assume that all the edges 

T, r', r", . . . are pointing to 7. As p G ker di, we have p^ + p^' + Pr" H = 0. 

Now as the cells a,a' ,a" , . . . ,a are ordered counter-clockwise around 7 and as 

the edges are pointing to 7, we have Sa,T = Sa'.r' = £a" .t" = • • • = 1, so that the 
sum 9i ([cr, a']) + 5^ ([cr', a"\ H V 3^ ([cr'", cr]) over a the boundary of a face 7* of 

r* is 0. 

By composition, for any loop of 7"*, the sum on the corresponding oriented edges 
is 0. This shows that the definition of q,, does not depend on the oriented path 
from cro to cr. 
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By construction, wc have 

For each interior edge t S T", there are two faces ai > (T2, which are adjacent to 
T. Thus, we have £0-1 ,t = —£(72,t and go-i = 90-2 + £<ti, tPt- We deduce that 

^ ^a^rQa^ ~ ^t7i,rQai ~t~ £(T2,t^(T2 ~ ^(T\,t{,Q(J2 ~t~ ^(TijtPt} ~t~ ^(J2,tQ(T2 ~ Pt' 

a 

This shows that 92(9) = P- In other words, im92 = ker^i and Hi{fRm,m' {T°)) = 
0. ' □ 

Proposition D.3. If Q has one connected component. Then -ff2(9^m,m'(T°)) = 

Proof. An element of H2{9{m,m' (T°)) = ker ^2 is a collection of polynomials {pa)aeT2 
such that G Rm.m' and p^r = Per' if f and a' share an (internal) edges. As T 
is a subdivision of a rectangle Dq, all faces a £ T2 share pairwise an edge. Thus 

Pa = Per' for all (7, a' e 72 and i?2(5^m,m' (T)) = Rm,m'- □ 

Notice that by counting the dimensions in the exact sequence $Hm,m'(7"), we 
recover the well-known Euler formula: /2 — /i + /o = 1 (the domain has one 
connected component). 



m 










; 






1 1 — 





















tp — i) 



(i — H'T 



4S f) 



(»— • 



— 



t3 




( 


1 










7 




h [ 


— 


i 




--E] 


h 












31 


( 


1 

'2 


S3 



^5 ( 
*4 ( 


1 •- 

1 i>- 

i -ii 






a 


( 

( 

i 


1 1> 

1 (1 

> E] 






1 ( •- 

1 •- 

Si Si 








( 

( 


1 ( 1 

1 <h 

is i'^ 


1 





— 


• ( 


1 


— • 


m 






(1 

^3 [3 

'^u 


1 ( 

1 1 


1 

7 
1 


— ( 1 

-f>--- 
-Hi 


H • 

E] 
-<( 






^, — 

Si 


• ( 

S2 i 


1 

3 


-• 

34 


-• 

S5 





